Self-localization of directed polymers and oppressive population control 
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We construct a phenomenological theory of self-localization of directed polymers in d+ 1 dimensions. 
In d = 1 we show that the polymer is always self-localized, whereas in d = 2 there is a phase transition 
between localized and free states. We also map this system to a model of population dynamics with 
fixed total population. Our previous results translate to static and expanding population clusters, 
depending on the birth and death rates. A novel "pseudo-travelling wave" is observed in some 
sectors of parameter space. 



PACS numbers: 05.40.-a, 87.23. Cc 

Directed polymers are important topological objects in 
many condensed matter systems. Examples are flux lines 
in superconductors |IJ , domain walls in ferromagnets , 
and atomic steps on vicinal surfaces Q . In the case of an 
isolated directed polymer, theoretical studies have con- 
centrated on the effects of quenched disorder, which may 
be either point-like or columnar. Using disorder to sup- 
press the wandering of flux lines is a central concept in 
designing useful high-temperature superconductors Q. 

Another mechanism by which the wandering of di- 
rected polymers may be suppressed is self-localization, 
mediated by elastic distortions of the medium in which 
the polymer is embedded. This has been studied previ- 
ously via a variational free energy functional which 
is the analogue of the quantum action for the polaron 
problem In the first part of this Letter we study 
the phenomenon from an alternative viewpoint, using the 
statistical mechanics of directed polymers. Such an ap- 
proach leads us to a simple phenomenological equation 
for the probability density P of the directed polymer. 
We show that in d = 1 the directed polymer is always 
self-localized, while in d = 2 the polymer undergoes a 
transition to self-localization above a critical value of the 
coupling constant A (which measures the strength of in- 
teraction between the polymer and the embedding envi- 
ronment). The details of this transition are investigated 
by numerical integration. 

In the latter half of this Letter we show that the phe- 
nomenological equation for P may be re-interpreted in 
the language of population dynamics. In this case P cor- 
responds to the local population density. The precise 
form of the dynamics is unusual, since along with the 
birth and death terms, we have a global constraint that 
fixes the total population. [Similar constraints have been 
studied recently in the context of plant population mod- 
els |?| , and a link has been made to self-organized critical- 
ity Also, interesting links have been made between "non- 
Hermitian quantum mechanics" and population biology 
H .] The issue in the present work is not whether the 
population thrives or becomes extinct, but rather how 
the individuals spatially distribute themselves. Different 
phases are observed corresponding to self-localized and 
expanding populations. We also test the predictions of 



the theory by simulations of a discrete reaction-diffusion 
model with a global constraint. Finally, within the con- 
text of population dynamics we are at liberty to reverse 
the sign of the coupling constant A. In this case, we find 
that in d = 1 the density is delocalized via a pseudo- 
travelling wave - an unusual scenario which interpolates 
between dynamical scaling and travelling waves. We end 
the Letter by noting the general mapping between di- 
rected polymers and constrained population dynamics. 

We consider a directed polymer with elastic constant 
k in a space with d transverse dimensions r, and one 
longitudinal dimension t. The polymer experiences ther- 
mal fluctuations from a reservoir at temperature T, and 
a potential V(r,t). One can write the partition function 
Z(r,t) for such a polymer with ends fixed at (0,0) and 
(r, t) in the form of a path integral. Following standard 
methods this path integral may be rewritten as a 
partial differential equation of the form 



d t Z 
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where v = T/2k and we have absorbed a factor of 
T into the potential. The probability density of the 
polymer is constructed by normalizing Z: P(r,t) — 
Z(r,t)/ J d d r' Z(r',t). Inserting this definition into the 
equation for Z yields a closed equation for P: 



d t P = vV 2 P - VP 



dVV(r',£)P(r',i) 



P. (2) 



To our knowledge, this description of the directed poly- 
mer has not been studied previously. Note this equation 
does not have the form of a Fokker-Planck equation - 
conservation of probability is not achieved via an equa- 
tion of continuity. The physics of directed polymers is 
intrinsically non-local. 

We want to consider the potential V to have been 
caused purely from elastic distortions of the embedding 
environment, due to the presence of the directed polymer 
itself. In principle we could formulate a detailed micro- 
scopic model which includes these elastic degrees of free- 
dom explicitly, and then attempt to integrate them out 
to find the effective potential. However, we are content 
to take a phenomenological viewpoint: the elastic distor- 
tions caused by the polymer will provide an attractive 
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potential which is strongest in the immediate vicinity of 
the polymer, and negligible far from the polymer. If we 
try to construct such a potential from the probability 
density P, we come at once upon the particularly sim- 
ple form V(r,t) = —XP(r,t), where A is a non-negative 
coupling constant. Inserting this relation into Eq.(^|) we 
have a closed equation for P of the form 

d t P = vV 2 P + XP 2 - Xa(t)P , (3) 

where a(t) = J d d r P{r,t) 2 . 

Given we have a non-linear integro-differential equa- 
tion for P there is little hope of a complete analytic solu- 
tion. However we will uncover much of the physics from 
scaling arguments, and numerical integration. First, let 
us consider the possibility of simple scaling, and make a 
similarity Ansatz, which, based on the radial symmetry 
of the problem, we take as P(r,t) — £(t)~ d F(z), with 
z = r/£. Note, the prefactor is fixed by normalization of 
P. Inserting this Ansatz into Eq.(||) yields 

i{dF + zF') + (v/0(F" + (d - l)F'/z) + 

(X/H d - 1 )(F-a )F = 0, (4) 

where (Jq = Sd J dzz d ~ 1 F(z) 2 , and Sd is the area of the 
unit hypersphere. We need to balance the terms arising 
from the derivative with respect to t and the Laplacian. 
This implies £ = (V/£) which in turn gives £ = [2vt) 1 / 2 . 
So for large t the correlation length £ is large, and we 
see on comparing the Laplacian terms to the interaction 
terms (i.e. those proportional to A) that for d > 2 the 
interaction terms are negligible. In other words, for d > 2 
the behaviour of P corresponds to thermal diffusion (at 
least for modest values of A) . 

With the assumption of a ^-dependent correlation 
length, we see that the interaction terms dominate for 
d < 2, and there is no self-consistent solution. We there- 
fore demand that £ is i-independent, and denote it as £o- 
Focusing on d — 1, we find 

F" + (X^M(F-a )F = . (5) 

This equation is simple enough to yield an exact solution. 
Expressed in terms of the original spatial variable x, we 
find the ^-independent solution 

P(x) = (l/4£o)sech 2 (z/2£ ) , (6) 

where the localization length is self-consistently found to 
be £o = 61//A. We have confirmed the stability of this 
solution by direct numerical integration of Eq. (||) using a 
variety of initial conditions. In all cases we find rapid con- 
vergence to the form given above. It is remarkable that 
our solution coincides exactly with the well-known form 
for the polaron probability density, first found by Rashba 
||. This agreement is satisfying from a physical view- 
point, but by no means obvious given the mathematical 



difference between Eq. (|3|) and the non-linear Schrodinger 
equation governing polaron physics. 

We can achieve a perfect balance of all the terms in 
Eq.([|) in exactly two dimensions. In this case, the cor- 
relation length drops out of the equation and we have an 
ordinary differential equation for F{z) of the form 

F" + (z + l/z)F' + 2F + X(F - <t )F = , (7) 

where A = Xjv is a dimensionless coupling constant. We 
have been unable to derive an analytic solution to this 
equation. We have studied its properties using a straight- 
forward numerical integration scheme, starting with a lo- 
calized initial condition, generally taken to be Gaussian. 
For modest values of A the probability density diffuses 
outwards, and we find good data collapse in accordance 
with the similarity Ansatz defined in terms of z = r / '\fi. 
As we increase A we find a sudden transition at A c ~ 30.5, 
above which the probability density shrinks to a micro- 
scopically collapsed state. Details of this transition will 
be given in a longer paper jl(J . 

We now move away from the physics of directed poly- 
mers, and re-interpret our central equation (0) within a 
very different context, namely reaction kinetics or popu- 
lation dynamics. Consider a species A (either chemical, 
biological or homo sapiens) which undergoes the follow- 
ing birth and death processes: A — > and 2A — > 3^4, with 
rates k and k' respectively, which are in general time de- 
pendent. We denote the relative density of A by p(x,t), 
and stress that the symbol t now represents real time. 
We may write a reaction-diffusion equation to describe 
these processes: 

d t p = DW 2 P + k'(t)p 2 -k{t)p , (8) 

where D is the effective diffusion constant of individual 
A 'particles'. The above equation is only an approximate 
description of this process, and a priori is only expected 
to be valid above some critical dimension d u . For lower 
dimensions, it is often necessary to explicitly include the 
appropriate (microscopically derived) noise term in the 
above equation [fll|| . For a qualitative understanding of 
the model (which is further refined below), such a noise 
term is not required. 

At this point, we shall introduce a somewhat unusual 
feature - a global constraint on the population. In other 
words we fix the normalization of the relative density 
to unity: j d d rp = 1. From Eq.(||) we see this im- 
poses k'{t)/k{t) — J d d rp 2 . We can imagine two special 
cases of the above. First, we can fix the death rate k 
to be constant, which implies that the total population 
is controlled by tuning the birth rate k' (which is not 
unknown). Alternatively, we have the chilling mecha- 
nism of total population control in which the birth rate 
is fixed, equal to ko say, and the death rate is adjusted 
according to k(t) = ko J d d xp 2 . At the level of our phe- 
nomenological description, these two choices lead to sim- 
ilar behaviour, so we shall concentrate on the second. 
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Substituting the appropriately tuned rates into Eq.(|^) 
we have 



dtp = DV 2 p + k oP 2 - k 



d d r'p{v',tf 



(9) 



This equation is identical in form to the self-consistent 
equation (j||) for the probability density of the directed 
polymer. Thus all of our results may be taken over and 
immediately applied to the population system outlined 
above. We state the gross features here. In d = 1 the 
population will distribute itself in a stationary localized 
structure, with a spatial scale £o ~ D/ko. In d = 2 the 
population will attempt to spread itself as widely as pos- 
sible, with a diffusive dynamics, but only for ko less than 
some critical value k c . Above this value, the system will 
undergo collapse and exist in a stationary high-density 
phase. For d — 3 the population will generally exist in 
the diffusive spreading phase (although there is a possi- 
bility for different behaviour for very high values of ko). 

As mentioned before, it is not obvious a priori whether 
a phenomenological description such as Eq.(^) is an ac- 
curate model of the true population dynamics for low di- 
mensions. Intuitively, we expect the microscopic fluctua- 
tions to be suppressed in our model due to the long-range 
interactions which are built-in via the global population 
constraint. In this case, our qualitative predictions from 
Eq.(Q) should hold true. We have tested this hypothesis 
by performing microscopic simulations of the birth/death 
process described above. Our model consists of N ran- 
dom walkers on a d-dimensional hypercubic lattice. The 
walkers obey exclusion statistics. If two walkers are at 
adjacent lattice sites, there is a probability q for them to 
spawn a new walker at a neighbouring unoccupied site. 
If this occurs, a walker is randomly selected from the en- 
tire population and is killed - thus conserving the overall 
population. We have studied this model in various di- 
mensions. In d = 1, and starting from a compact clus- 
ter of 100 particles, we find that the cluster spreads out 
somewhat, but that the mean cluster size asymptotes to 
a constant. On averaging over realizations we find that 
the density profile of the cluster (in the center of mass 
frame) is in approximate agreement with the solution of 
the continuum model. We refer the reader to Fig.l where 
a direct comparison is made. 

In d = 2 we find a phase transition on varying the 
spawning probability q. For q < q c the cluster asymp- 
totically diffuses outwards, with a cluster size growing as 
y/i. For q > q c the cluster remains compact. The value 
of q c is generally very small. For example, for N = 400, 
we find q c ~ 0.004. The precise N dependence of q c , 
along with details of the simulations will be given in a 
longer work |u|. At a qualitative level, we see that the 
continuum formulation (|9|) provides a good description 
of the underlying microscopic model. 




FIG. 1. Average (stationary) cluster density from simula- 
tions (solid line) in d — 1 with N = 100 particles, compared 
to the continuum prediction (dashed line) with £o fitted 
to the mean cluster size. 

In the directed polymer context there is no compelling 
reason to study the system for A < 0. However, within 
the context of population dynamics, this is a perfectly 
acceptable sector of parameter space to explore. In order 
to motivate such a model, we shall consider a different 
birth/death process to the one described above. Con- 
sider now that the individuals A interact with the rules: 
A — ► 2 A and 2 A — > 0, with rates r and r' respectively, 
which are in general time dependent. In other words, 
we have rules corresponding to mutual annihilation and 
single-sex reproduction. We choose to take r' — r (con- 
stant) and tune r to maintain a constant total popula- 
tion. The phenomenological reaction-diffusion equation 
for this process takes the form 



d t p = £>V 2 / 



rop 



ro 



d d r'p(r\tf 



P 



(10) 



As desired, this corresponds to our original directed poly- 
mer equation (|3j) but with a negative coupling constant. 

This equation has some remarkable properties, which 
may be uncovered from an approximate analytic treat- 
ment. We shall restrict our attention to d = 1. We 
note that spatial coupling enters the model in two dis- 
tinct ways: locally via the Laplacian, and globally via 
the normalization constraint. This means that there can 
be non-trivial dynamics even if we set D — - a limiting 
case which allows an exact solution of the equation. We 
shall confine ourselves to stating a few pertinent results. 
Details will be given in a longer work p0| . It is obvious 
that in this limit the dynamics will be extremely sensitive 
to the initial condition. One static solution corresponds 
to choosing the initial condition to be a normalized top- 
hat function. More interesting solutions are obtained 
with initial conditions with infinite support. Consider 
p(x,0) = (j4/£o) ex p(— (M/£o)^)- For long times, one 
finds that this profile spreads outwards, with a flat cen- 
tral region of height ~ l/£(i), where is the lateral 
size of the profile. We find that this length scale grows 
in time as £(i) ~ t a with a = 1/(1-1-/3). The solution for 
the density profile is expressed as 

p( X ,t)^(i/m)F((\x\/^r-m)/^r) , (n) 
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where F(y) is simply related to the initial profile. Since 
this solution does not have a canonical dynamical scal- 
ing form, there may be other relevant scales; and indeed 
we find that there is a second important length scale, 
which is the width W(t) of the interface separating the 
region of "flat" nonzero density and the region of van- 
ishing density. This quantity changes in time accord- 
ing to W(t) ~ COXVCW)' 3 - Therefore the interface 
width shrinks (grows) with time for f3 > 1 (/? < 1). The 
above solution is rather unusual as it, in some way, in- 
terpolates between a travelling wave solution F(x — vt) 
and a scaling (or similarity) solution F(x/£(t)) We 
term it a "pseudo-travelling wave". [We note in pass- 
ing that an initial condition with tails decaying with a 
power q evolves under normal dynamical scaling with 
p ~ (l/t)F(\x\/t) and F(z) ~ for z > 1.] 



0.001 



0.0008 



0.0006 



0.0004 



0.0002 




numerical integration of Eq. ( |l0| ) in d = 1 for a variety of 
initial conditions, and with zero or non-zero D. The ex- 
act results for D = have been numerically confirmed, 
however the situation for D > is less clear. In Fig. 2 
we show an example of the evolving density grown from 
an initial Gaussian profile, with D = 1 and ro = 5. In 
the inset we plot the profile size and interfacial width. 
These quantities exhibit very strong corrections to sim- 
ple power law scaling, but a closer analysis (studying the 
'running exponents') indicates that the asymptotic power 
laws are approximately 0.65(1) and 0.30(1) respectively, 
which are close to the values obtained above. 

In conclusion we have considered two very different 
physical scenarios which share a common mathematical 
description. These are self-localization of a directed poly- 
mer, and birth/death models with fixed population num- 
ber. The physics of these processes is non-local which 
leads to a number of interesting results, including phase 
transitions from localized to non- localized states, and 
pseudo-travelling waves. 

Finally, it is noteworthy that the fundamental directed 
polymer equation (|^) with a general external potential 
V can be recast in terms of constrained population dy- 
namics. The case of a random potential is especially in- 
teresting, and corresponds to local, random birth/death 
processes, along with a global feedback to control the 
population. We are currently using this mapping to gain 
new insights into these two fascinating problems. 




FIG. 2. Density profiles for three different times, from nu- 
merical integration of ( jlo| ) with D = 1 and ro = 5. The inset 
shows a log-log plot of the time evolution of the cluster size 
£ and the interface width W . The straight lines (which are a 
guide to the eye) have slopes of 2/3 and 1/3 respectively. 

One may now ask how the above picture is modified 
when one re-instates the local coupling D. We have sev- 
eral arguments which all agree that, so long as the initial 
profile has a finite transverse scale, then for late times 
the density will have the form (|ll|), but with a selected 
value of f3 — 1/2. This implies that the profile grows out- 
wards with a transverse scale £(t) ~ t 2 ^ 3 and interfacial 
width W(t) ~ t 1 / 3 . The simplest (albeit non-rigorous) 
way to see why this value of (3 is selected is to assume 
that for non-zero D the density profile has a form approx- 
imately given by (|ll|), and to then ask how the Lapla- 
cian term in ( fj"o| ) can balance with the reaction terms. 
Since the profile is flat for |x| <C and vanishing for 
|x| ^> £(i) it is clear that the size of the Laplacian is 
set by the width of the interfacial region. We therefore 
estimate d 2 p ~ Similarly the reaction terms 

are of order p 2 ~ l/£ 2 - Balancing the two leads us to 
W{t) ~ £(i) 1/ ' 2 , which on comparing with our previous 
result W(t) ~ £(i) 1_/3 gives /3 — 1/2. We have performed 
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